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ABSTRACT 

The inner regions of accretion disks of weakly magnetized neutron stars are affected 
by general relativistic gravity and stellar magnetic fields. Even for field strengths 
sufficiently small so that there is no well-defined magnetosphere surrounding the 
neutron star, there is still a region in the disk where magnetic field stress plays an 
important dynamical role. We construct magnetic slim disk models appropriate for 
neutron stars in low-mass X-ray binaries (LMXBs) which incorporate the effects of 
both magnetic fields and general relativity (GR). The magnetic field — disk interaction 
is treated in a phenomenological manner, allowing for both closed and open field 
configurations. We show that even for surface magnetic fields as weak as 10'' — 10^ G, 
the sonic point of the accretion flow can be significantly modified from the pure GR 
value (near tgr = 6GM/c^ for slowly-rotating neutron stars). We derive an analytical 
expression for the sonic radius in the limit of small disk viscosity and pressure. We 
show that the sonic radius mainly depends on the stellar surface field strength Bq 
and mass accretion rate M through the ratio 6^ oc (3Bq/M, where /? ~ \B(p/Bz\ 
measures the azimuthal pitch angle of the magnetic field threading the disk. The 
sonic radius thus obtained approaches the usual Alfven radius for high b'^ (for which 
a genuine magnetosphere is expected to form), and asymptotes to QGM/c^ as iP' 0. 
We therefore suggest that for neutron stars in LMXBs, the distinction between the 
disk sonic radius and the magnetosphere radius may not exist; there is only one 
"generalized" sonic radius which is determined by both the GR effect and the magnetic 
effect. 

We apply our theoretical results to the kHz quasi-periodic oscillations (QPOs) 
observed in the X-ray fiuxes of LMXBs. If these QPOs are associated with the 
orbital frequency at the inner radius of the disk, then the QPO frequencies and their 
correlation with mass accretion rate can provide useful diagnostics on the (highly 
uncertain) nature of the magnetic field ~ disk interactions. In particular, a tight upper 
limit to the surface magnetic field Bq can be obtained, i.e., ^ 3 x 10'^(Mi7//3)i/2 Q, 
where Mn = M/(10^^ gs~^), in order to produce kHz orbital frequency at the sonic 
radius. Current observational data may suggest that the magnetic fields in LMXBs 
have complex topology. 
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1. Introduction 

The inner region of disk accretion onto neutron stars may be characterized by two unique 
radii: (i) The marginally stable orbit due to general gravity (GR). For nonrotating neutron stars 
this is located at QGM 

rcR = = 12.4A/1.4 km, (1) 

where M is the neutron star mass, and Mi. 4 = M/^IAMq). For finite rotation rates, tgr is 
somewhat smaller. The fiow behavior near tgr has been subjected to numerous studies, especially 
in the context of black hole accretion disks (e.g., Muchotrzeb &; Paczynski 1982; Matsumoto 
et al. 1984; Abramowicz et al. 1988; Narayan et al. 1997; Chen et al. 1997): Close to roR the 
inward radial velocity of the accreting gas increases steeply with decreasing radius and becomes 
supersonic. The existence of such marginally stable orbit for neutron star is predicated on the 
fact that neutron star models constrvictcd using different nuclear equations of state generally 
give a stellar radius less than tgr (Arnett & Bowers 1977; Kluzniak &: Wagoner 1985). (ii) The 
magnetospheric radius, rm, below which magnetic stress dominates disk plasma stress. While 
the precise value of depends on the (rather uncertain) details of the magnetic field - disk 
interactions, it is estimated to close to or slightly less than (by a factor of a few) the spherical 
Alfven radius, i.e., 

mJgMr ) = ^^vR^M^l^Bfu^?' (km), (2) 

with r? <, 1 (e.g., Pringle & Rees 1972; Lamb, Pethick & Pines 1973; Ghosh & Lamb 1979; 
Arons 1987), where we have scaled various quantities to values appropriate for neutron stars in 
low-mass X-ray binaries (LMXBs): R = 10i?io km is the neutron star radius, Bq = 10^ G 
is the dipolar surface field strength, and M = (10^''g s~^)Mn is the mass accretion rate (The 
Eddington accretion rate is about 10^^ g s~^). For highly magnetized neutron stars (such as X-ray 
pulsars, typically having B ^ 10^^ G), is much greater than tgr and the stellar radius, the 
disk is therefore truncated near rm, within which the disk plasma becomes tied to the closed field 
lines and is funneled onto the magnetic poles of the star, although some plasma may continue to 
fall in the equatorial plane as a result of interchange instabilities (Spruit &; Taam 1990; see also 
Arons & Lea 1980). For weakly magnetized neutron stars, such as those expected in LMXBs, v-^ 
and rGR are comparable, and the plasma may not climb onto the field lines before reaching the 
stellar surface. A question therefore arises as to how the magnetic field affects the the dynamics 
of the inner disk and changes the sonic point. In this paper, we present an unified (albeit 
phenomenological) treatment of neutron star accretion disks under the combined infiuences of 
magnetic fields and strong gravity. 
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Our study is motivated by the recent observations using the Rossi X-ray Timing Explorer 
(RXTE) (Bradt, Rothschild & Swank 1993) which revealed kilo-Hertz quasi-periodic oscillations 
(QPOs) in the X-ray fluxes of at least thirteen LMXBs (see Van dcr Klis 1997 for a review). These 
kHz QPOs are characterized by their high levels of coherence (with v j up to 100), large rms 
amplitudes (up to 20%), and wide span of frequencies (500 — 1200 Hz) which, in most cases, are 
strongly correlated with the X-ray fluxes. In several sources, the X-ray power spectra show twin 
kHz peaks moving up and down in frequency together, with the separation frequency roughly 
constant. Moreover, in five atoll sources single QPOs (with a much higher level of coherence) have 
been seen during one or more X-ray bursts, with frequencies equal to the frequency differences 
between the two peaks or twice that. This is a strong indication of beat phenomena (Strohmayer 
et al. 1996). While the origin of these QPOs is uncertain, it is clear that the action must take 
place close to the neutron star, either in the accreting atmosphere (Klein et al. 1996) or in the 
inner disk (Strohmayer et al. 1996; Miller, Lamb and Psaltis 1996). A generic beat-frequency 
model assumes that the QPO with the higher frequency is associated with the Kepler motion at 
some preferred orbital radius around the neutron star, while the lower-frequency QPO results from 
the beat between the Kepler frequency and the neutron star spin frequency. It has been suggested 
that this preferred radius is the magnetosphere radius (Strohmayer et al. 1996) or the sonic radius 
of the disk accretion flow (Miller et al. 1996). 

In this paper, we are not concerned with the actual mechanisms by which kHz QPOs in 
the X-ray fluxes of LMXBs may be produced (See Miller et al. 1996 and Van der Klis 1997 for 
extensive discussion on various possibilities). Rather, our main purpose is to understand what 
physical effects determine the characteristics of the inner accretion disks in LMXBs. In the 
sonic-point model. Miller et al. (1996) suggest that some accreting gas can penetrate inside the 
magnetosphere, whose boundary is located at a larger radius than the sonic radius. For unknown 
reasons, they assume that these gases are unaffected by the magnetic field once they are inside the 
magnetosphere and remains in a Keplerian disk. They further suggest that the variation of QPO 
frequency results from the change in radiative forces on the accretion disk. We note, however, that 
the effect of radiative forces on the disk fluid (rather than test particle orbiting the central star), 
is far from clear. Calculating particle trajectories without solving for the global disk structure 
(M. C. Miller 1997, private communication) is inadequate for determining the magnitude of the 
radiative forces. While the radiative forces may be important for high-luminosity Z-sources, their 
effects on the disk dynamics are expected to be be small for low-luminosity systems (L less than 
10% of i^Edd)- On the other hand, it is well known that millisecond pulsars have magnetic fields 
in the range of 10^ — 10^ G, and one expect that neutron stars in LMXBs to have the similar 
range of field strengths. While the magnetic field may not be strong enough to induce a corotating 
magnetosphere outside the neutron star, it can nevertheless influence the dynamics of the inner 
disk flow by transporting away angular momentum from the disk. 

Ideally, to properly assess the dynamical effect of magnetic fields on the accretion disks, one 
needs to solve for both the fluids and the fields self-consistently. This is a difficult task if not 
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impossible. Despite many decades of theoretical studies (e.g., Pringle &: Rees 1972; Lamb et 
al. 1973; Ghosh & Lamb 1979; Aly 1985, 1991; Arons 1987; Spruit & Taam 1990; Sturrock 1991; 
Shu et al. 1994; Lovelace et al. 1987, 1995; Stone & Norman 1994; Miller & Stone 1997), there 
remain considerable uncertainties on the nature of the stellar magnetic field - disk interactions. 
Particularly outstanding are the issues related to the efficiency of magnetic field dissipation in 
and outside the disk and whether the stellar field threads the disk in a closed configuration or it 
becomes open due to differential shearing between the star and the disk. It seems unlikely that 
some of these issues can be resolved on purely theoretical grounds. In this paper, we shall not 
attempt a self-consistent magnetohydrodynamics (MHD) calculations. Rather, we shall adopt a 
phenomenological approach and consider rather general field configurations. We believe that such 
an approach is useful in bridging the gap between full MHD theories and observations. Indeed, 
as we show in this paper, if the observed kHz QPOs are associated with the sonic point Kepler 
frequency, then various systematics of kHz QPOs should provide useful constraints on the nature 
of magnetic field - disk interactions as well as on the magnetic field structure in LMXBs. 

In §2 we introduce a model of magnetic slim accretion disk. Numerical solutions are presented 

in §3. Because of various uncertainties in disk parameters, we shall focus on the simplest models, 
leaving more complete exploration to future studies. However, in §4 we introduce the notion of 
"generalized marginally stable orbit" including both the GR and magnetic effects. We derive an 
analytical expression for the sonic radius, which shows that the sonic point depends mainly on two 
parameters characterizing the disk magnetic field (in addition to the neutron star mass). We show 
how different modes of magnetic field - disk interactions can lead to different sonic-point orbital 
frequencies and their scalings with the field strength and mass accretion rate. Section 5 concerns 
the equilibrium spin periods of neutron stars in our slim magnetic disk model. Some applications 
to kHz QPOs in LMXBs are discussed in §6, where we show that our phenomenological approach 
can be used to learn about the physics of magnetic field - disk interactions. 

Unless otherwise noted, we use geometrized units in which the speed of light and Newton's 
gravitation constant are unity. 



2. Slim Magnetic Accretion Disks 

We now consider geometrically thin axisymmetric accretion disk in steady state, taking into 
account of the transonic nature of the flow, and the deviation from Keplerian motion in the 
inner region of the disk. Our models generalize the usual "slim disks" around black holes (e.g., 

Muchotrzeb & Paczyhski 1982; Matsumoto et al. 1984; Abramowicz et al. 1988; Narayan et 
al. 1997; Chen et al. 1997) by including the effect of magnetic fields. GR effects are included in our 
purely Newtonian treatment by using the pseudo-Newtonian potential introduced by Paczyhski &: 
Wiita (1980) 

M 
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This potential correctly reproduces the marginally stable orbit Q located at tgr = 6M, and is 
adequate for this initial exploration, considering the much greater uncertainties in the magnetic 
field - disk interactions. The self- gravity of the disk is neglected. 

We assume that the accreting material is confined to a thin disk, and we do not formally 
introduce a magnetosphere in our model. As discussed in §1, when the field strength is sufficiently 
high (as in X-ray pulsars, which typically have B ^ 10^^ G), there is no question that a corotating 
magnetosphere exists outside the neutron star surface, located near (see Eq. note that, 
theoretically, the magnetosphere radius is not known to within a factor of two, nor is it clear 
what the width of the transition zone is). In such a high-field regime, we shall find that the sonic 
point as obtained from our model is approximately equal to the usual Alfven radius. Although in 
our model the fiow continues to be confined in the disk plane even inside the magnetosphere, in 
reality it may well behave differently (e.g., the plasma may "jump" onto the field lines and get 
funneled onto the magnetic poles). Thus for high magnetic systems, the supersonic portion of our 
fiow (inside the sonic point) may not be realistic. However, for low magnetic systems (such as 
LMXBs), which is the main focus of this paper, there needs not be a genuine magnetosphere to 
truncate the disk flow, but the magnetic forces can still shift the sonic point to a radius larger 
than 6M. In such low-fleld regimes, we expect our global flow solutions to have a wider validity. 



2.1. Basic Equations 

The mass continuity equation takes the form 

M = -27rrE'u, (4) 

where u is the radial velocity of the flow (n < for accretion), and S = Jdz p ~ 2Hp is the surface 
density of the disk. The disk half-thickness is given hy H = Cg/i^K, where Cs = (p/p)^^'^ is the 
isothermal sound speed, and is the Keplerian angular velocity (for the pseudo-Newtonian 
potential): 

The radial momentum equation reads 

W P _ _ R R 

(6) 



dr dr ' 27rS 



^With small spin parameter a = J/M^ << 1 (where J is the spin angular momentum), we have roa/M ~ 
6(1 - 0.544a), and AmK{rGR) ^ 6~^^^{l + 0.748a). For a spin frequency of 300 Hz (Strohmayer at al. 1996), this 
amounts to a correction of 7% to tgr and 10% to These corrections are neglected in this paper. 

^We neglect the magnetic field effect on the disk thickness. In reality, the disk can be compressed or flattened 
depending upon the field configuration (see Wang et al. 1990 and Stone & Norman 1994). However, the subsequent 
analyses in this paper do not depend upon having an explicit expression for H since only height-integrated equations 
are used. 
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where P = Jdzp is the integrated disk pressure, O is the angular velocity. The last term in Eq. (|6|) 
represents the dominant radial magnetic force, obtained by integrating over height the force per 
unit volume d{BzBr/^n)/dz and dividing by S. Note that in Eq. (|6|), Br\z=H is evaluated at the 
upper disk plane, and Br\z=-H = —Br\z=H- In deriving the magnetic force, we have also neglected 
the rr-component of the magnetic stress. 

The angular momentum equation reads 

dl r 



1 d / 3^ B,B^ 



.r^ dr \ dr J ' 2tt z=h\ ' ^ ^ 

where the second term on the right-hand-side is the magnetic torque per unit mass, obtained by 
integrating over height the torque per unit volume rBz{dB^/ dz) / At: and dividing by S. Equation 
can be integrated in r to give the conservation equation for angular momentum: 

MlQ = Ml + 2TTur^T.'^ + MNB{r), (8) 
dr 

where Iq is the integration constant, and 



/■oo 

MNsir) = - drr'^BzB^ 

Jr 



(9) 

The three terms on the right-hand-side of Eq. (^) correspond to advective angular momentum 
transport, viscous torque, and magnetic torque due the threading field lines from r to oo, 
respectively. The constant Iq is the eigenvalue of the problem; it should be determined by requiring 
the flow to be regular at the sonic point. We shall adopt the standard a— prescription for the disk 
kinematic viscosity, i.e., v = aHcg (Shakura & Sunyaev 1973). 

The energy equation describing the thermal state of the flow can be written in the form: 

^Tu— = ^visc + ^Joule - 2Fz. (10) 

dr 

Here S is the specific entropy (per unit mass), -Evisc is the viscous heating rate per unit area. With 
the a-prescription for the disk viscosity, we have 

1 / dn\'^ ( dVL\^ , , 

= (p..-) =.E(r-) . (11) 

The vertical (optically thick) radiative transport flux is 

ipK dZ KZj 

where k is the opacity and aT^ is the radiation energy density. The Joule heating rate -Ejouie 
depends on the field dissipation in the disk, and its specific form depends on our ansatz for the 
magnetic field (§2.2). 

Finally we need equations of state. For the inner disk region of interest in this paper, radiation 
pressure dominates over gas pressure. Thus we have p = and P = 2HaT^/?,. Also, the 

opacity is dominated by Thomson scattering, n = 0.4 cm^/g. 
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2.2. Ansatz for the Magnetic Field 

The equations above can be apphed to general axisymmetric magnetic field - disk 
configurations, as long as mass loss from possible disk wind is negligible, and the accreting material 
is confined to a thin disk plane. We now specify our ansatz for the magnetic fields. The vertical 
field component is assumed to take the form 

B, = Boi-] . (13) 



r 

We shall mostly focus on the n = 3 case, corresponding to a central stellar dipole field threading 
the disk (e.g., Ghosh & Lamb 1979; Konigl 1991; Yi 1995; Wang 1995), ahhough we will also 
consider more general values of n, as in the cases when high-order multipoles are important (Arons 
1993) or when field lines become open due differential shearing between the disk and the star (Aly 
1985, 1991; Sturrock 1991; Newman et al. 1992; Lynden-Beh & Boily 1994; Lovelace et al. 1995). 
In reality, the power-law relation in Eq. ([l^ ) is most likely to be valid only for a small range of 
r, but as we shall see in §4, the sonic point is mainly determined by the local behavior of the 
magnetic field. 

For the azimuthal component of the magnetic field, we consider two possibilities: 

(i) If the stellar magnetic field threads the accretion disk in a closed configurations (e.g., 
Ghosh & Lamb 1979), then B^j, is governed by dB^/dt = Bzd{rQ)/dz — B^/t (where r is the field 
dissipation time). In steady-state, this gives B^ ~ t{^}s — / H)Bz, where is the rotation 
frequency of the star. We define a dimensionless parameter (3 such that 

B, =P (%^) S., (14) 

z=H \ ^^K ) 

where is the Kepler frequency and r2(~ ^k) is the disk orbital frequency. Various (uncertain) 
estimates for the field dissipation timescale in the magnetically threaded disk configurations have 
been summarized in Wang (1995). 

(ii) If the magnetic field becomes open (e.g., Lovelace et al. 1995), we assume 

B^\z=H = -PB„ (15) 

where /5 ~ 1 specifies the maximum twist angle of any field line connecting the star and the disk. 
Clearly, Eq. ( [T^ ) encompasses the second possibility if we set 0^ = 0. However, we note that the 
physical meaning of Iq in these two cases are rather different: For closed field configurations (i), 
MIq measures the total torque on the star, while for the (partially) open field configurations (ii), 
MIq also include the angular momentum carried away from the disk by the magnetic fields of disk 
outflow. In both cases, MIq is the total angular momentum transported away from the disk per 
unit time. 



Similar to Eq. (|l^, our ansatz for the radial component B^ of the disk magnetic field is 

(— 

z=H \^Kr 



Br 



Pr(TT^] Bz. (16) 
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We expect Pr to be of the same order of magnitude as /3. 

With the particular ansatz for the magnetic fields given by Eq. (p^, the Joule heating rate 
-f'jouie can be calculated as 

= /d. ^^B.B,§-m = i^PrBi^-^. (17) 
The dissipation due to the current associated with Br is much smaller and has been neglected. 



3. Numerical Solutions 

It is useful to define dimensionless field strengths b and br via 

= = 7.317 PRfM-r^, (18) 

b'r = Pr^ = 7.317PrR%'M[r^, (19) 

where Ir = {MR)^/'^ . Roughly speaking, 5^ is the ratio of the total magnetic torque and the 
characteristic accretion torque on the neutron star. Comparing with Eq. (||) we see that for dipole 
magnetic fields, ta/R ~ {b'^/Pf^^ ~ b^^^ In all our calculations, we choose br = b; The flow 
structure and the sonic radius are rather insensitive to the value of br- The radial force equation 
and the continuity equation can be cast in the form which reveals the existence of a sonic point: 

where = (dP/(iS)fiow The sonic point (where |m| = a^) is a critical point of the differential 
equation. The other equations can also be rewritten in the forms convenient for numerical 
integration. The angular momentum equation is 



with 

dr r \r J 

The energy equation is 



M dr IttuVLk \dr J 2tt r'^flx Kr J kM 



The eigenvalue Iq is adjusted so that the solution is regular at the sonic point. Equations 
(p0|)-(23) are integrated inward from an outer radius (far from the sonic point) where the disk is 
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approximately Keplerian. Note that in this outer Keplerian region, Eq. (22) can be integrated to 
give 



NB{r) 



b'^lR /i?^2n-3 



1 



(4n - 6)17, 



1 



(8n - 18)M 
(4n — 7)r 



(24) 



2n-3\rJ [' (4n - 9)(M/r3)i/2 
and the angular momentum equation yields 

where Ik = r'^^K- The sound speed Cg can be obtained from the energy equation (|23|). Since the 
radial velocity is small at large radius, the entropy advection term can be neglected. Substituting 
ac^/u from Eq. (^) into Eq. (^), we find 



kMQk 
Attc 



(26) 



Equations (^OD-(^) turn out to be a rather stiff set of equations. We have opted to adopt 
a further simplification by assuming the disk is isothermal. We estimate the range of using 
the thin disk expression ( |26| ) evaluated near the sonic radius. The physical rationale behind this 
approximation is that the transition from Keplerian disk to supersonic flow is very sharp, and 
we do not expect the sound speed to change significantly in this transition region (near the sonic 
point). As we shall see in §4, the sonic radius is insensitive to the thermal state of the disk 
when the sound speed is small. We have not studied the general dependence of the sonic radius 
on the thermal state of the disk. However, considering the very large uncertainties in the disk 
magnetic fields, the isothermal approximation should be adequate for use in the first step in our 
investigation. 

Figure 1 depicts two examples of transonic accretion fiows, with closed dipolar stellar fields 
threading the disks (n = 3). We choose a standard set of parameters: a = 0.1, Cg = 0.01, 
R/M = 5 (a typical value for a M = 1.4M0, i? = 10 km neutron star), and = 0.013M-1 
(corresponding to stellar spin frequency of 300 M-,"^ Hz). As expected, the magnetic fields slow 
down the tangential flow velocity, and, together with the strong relativistic gravity, make the 
radial velocity supersonic at small radii. For small 6^, the deviation of O from Qk is small; for 
larger 6^, the sonic radius larger, and Q gradually approaches 17, inside the sonic point. In 
these examples, the sonic points are located at r,/M = 7.43 (for = 1) and 10.95 (for = 10), 
the corresponding eigenvalues (/q) are 3.850M and 4.153M, respectively. 

In Figure 2 we show the sonic radius and the corresponding speciflc angular momentum Iq 
as a function of b^ for several different values of a and c,. We assume R/M = 5 and Qs = 0.013/M 
for these models. The following trends have been found: For a given 6^, a larger a tends to make 
r<j larger (i.e., viscosity tends to "destablize" the disk), while a larger Cg tends to make r<j smaller 
(i.e., pressure "stablizes" the disk). However, we emphasize that dependences of on these disk 
parameters (a,c<j) are rather weak. Moreover, as Cg decreases, the sonic converges to a value 
independent of a and Cg — The reason for this convergence will become clear in §4. 
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4. Analytical Approximation to the Sonic Radius 

As our numerical results in §3 indicate, in the limit of small disk pressure and viscosity, the 
sonic radius approaches a value independent of the disk parameters (a, Cg and the equation of 
state). This asymptotic sonic radius can be derived analytically using a simple mechanical model: 
Consider the equation of motion of a test mass around a neutron star 

|^ = ^-f^ir, k = l + NB{r). (27) 

Imagine that the test mass is "attached" to a magnetic field line so that its orbital angular 
momentum / is not conserved by itself. The conserved angular momentum Iq includes the 
contributions from both the orbital motion and the angular momentum A'^^ [cf. Eq. ^] carried by 
the threading magnetic field. Equation (|27|) can be obtained by setting the pressure and viscosity 
to zero in our general slim disk equations (§2.1). The radial component of the magnetic force has 
been neglected. An equilibrium orbit is determined by the condition 

Io-Nb = n\r = lK{r). (28) 

Deviation Ar from the equilibrium is governed by the perturbation equation of the form 

^ + ^2^Ar = 0, (29) 
where the the effective epicyclic frequency Keff is given by 

(Recall that In = y/ MR and R is the neutron star radius). Setting Kefj = 0, we obtain a critical 
orbit, which we dub the "generalized marginal stable orbit", located at r = tmso- Clearly, tmso 
depends only on the gravitational potential and the local magnetic torque (dNB/dr) = r'^BzB^/M. 
The magnetic field enters only through the dimensionless ratio b. The corresponding constant 
eigenvalue = Ik (rMSo) + -^b(^mso)) however, depends on the global field structure. 

In Figure 2 we plot tmso and Iq against 6^ for = 0.013/M (corresponding to spin frequency 
of 300 M{;] Hz) and R/M = 5. We see that tmso is the upper limit to the numerically determined 
r<j, i.e., tmso the asymptotic sonic radius as the disk viscosity and pressure diminish. 

It is of interest to consider two limiting cases: (i) In the Newtonian limit (neglecting the GR 
effect), or equivalently when 6^ is large (so that tmso >> 6M), we have 



^MSO 

R 



2bHl ^ 



2/(4n-5) 

(31) 



For n = 3, this is the standard result for the inner radius of the Keplerian disk, as determined by 
M{dlK/dr) = -r'^B^B^l^^jj (e.g., Arons 1993; Wang 1995); (ii) In the limit of smaU 6^ (so that 
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'"MSO is close to 6M), we find 



M 



16 &2 /i^x(4n-5)/2 r 





(32) 



6 + 



This gives the correction to the standard general-relativity-induced MSO located at tgr = 6M. 

The consideration of the limiting cases clearly indicates the sonic point (or the generalized 
MSO) constructed in our slim disk model includes the essential physics embodied in the 
determination of the usual magnetosphere radius. As we argued at the beginning of §2, for high 
magnetic systems, a genuine magnetosphere should certainly exist outside the neutron star. For 
low magnetic systems, however, the accretion flow may well be conflned in the disk, and the 
distinction between the sonic point and the magnetosphere boundary may not exist. 

We emphasize that our analysis given here is phenomenological. It only takes account of the 
dynamics of disk under a fixed magnetic field configuration, while a full MHD treatment should 
include perturbations of both disk fiuid and magnetic fields. The usefulness of our analytical 
result lies in the fact that tmso provides a good approximation to the sonic radius; and the sonic 
point is induced by both the GR effect and the magnetic effect. We note that in the presence of 
disk viscosity and magnetic fields, a fluid element continuously falls inward by the viscous stress 
and the magnetic stress, and therefore the concept of "marginally stable orbit" does not strictly 
apply. Nevertheless, we use the term "generalized MSO" to refer to the asymptotic sonic radius 
as determined by our analytical expressions. 

Figure 3 depicts the orbital frequency ^ at the generalized MSO (approximately the sonic 
radius) as a function of the mass accretion rate M for several different combinations of model 
parameters (n and fi^). Note that once we specify n and in units of 1/M, the numerical value 
of ry[so/M depends on the other parameters only through the combination 6^(i?/M)'^^"~^)/^ (see 



as the x-variable in Fig. 3 (where By is the surface field Bq in units of 10^ G). For large M /{(3Bj) 
(or small 6^), tmso approaches tgr = 6M and j^mso = '^ii'('^MSo) approaches 1.57/Mi,4 kHz. The 



shall discuss some of the applications of Fig. 3 to QPOs in LMXBs in §6. 



^To be consistent with the relativistic expression, we use umbo = (M/rMgo)^'^^/(27r), rather than the approximate 
pseudo-Newtonian Eq. (|5|). 



Eqs. We have therefore used 




(33) 



scaling of z^mso with M depends mainly on n, the index which specifies the shape of the magnetic 
field lines (see Eq. [Q). Note that the sonic point converges to the MSO only in the limit of 
small Cs- Thus in reality, the dependence of the sonic-point Kepler frequency i'Ki'>'s) on M may 
be different if Cg is not small. For example, Cg (x M for radiation-dominated optically-thick disk, 
and therefore the dependence of vxifs) on M is slightly steeper than what is shown in Fig. 3. We 
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5. Equilibrium Spin Period of Neutron Star 



As discussed before (§2.2), for closed field configurations, the quantity MIq measures the total 
torque on the neutron star due to the accreting matter and the threading magnetic fields. It is of 
interest to consider how the equilibrium stellar rotation rate ^s,eq, at which /q = 0, is determined 
in the slim disk model. We shall restrict to the dipole fields (n = 3) in this section. 

First consider the result when the GR effect is neglected. In this case the Keplerian 



disk boundary is determined by the condition M^dlx/dr) = 
MlK{rm) = 2rl^(3Bl{rjn){l - Ws), where lOs = Vts/VtK{rm)- We find 



\z=H-' 



or 



2 62(1 



(cf. Eq. pT|). The total torque on the star is given by 



iVt,, 



ML 



MlK{rm)l 





2/7 



1 - (8/7)a;, 



(34) 



(35) 



Equilibrium is obtained for LOg = 7/8. Using Eq. (p3) we then find = {b'^/Afl'^R and 



s,eq 



(7/8)(M/i23)i/2(52/4)-3/7^ ^Y,ich 



gives 



3.44Mi'//i?io 



-3/7 



3/2^-6/7 (j^jjz) = 1.467 Mf//i27o^®/^Mf/^ (/^^l) ' (kHz) 



(36) 



This is the standard result that the equilibrium spin frequency is equal to the Keplerian frequency 
at the Alfven radius. This result is plotted in Fig. 4. 

In the slim magnetic disk model, the total torque on the neutron star MIq is obtained from the 
eigenvalue Iq. In the asymptotic regime discussed in §4, can be determined from the analytical 
expressions: is obtained from the condition = in Eq. (|^), and then Iq = Ixifs) + A^s(rs) 
with Nb given by Eq. ( p^ ) (specialized to n = 3). It is straightforward to show that in the limit 
of r >> M, we recover the results given in Eqs. (34-^6|). The equilibrium ^ls,eq is obtained by 
requiring Iq = 0. In Fig. 4 we plot the equilibrium spin frequency function of 6^. Clearly, 

for large iP', our calculation agrees with the usual nonrelativistic result. Relativistic corrections 
are significant only when 6^ is small, for which the sonic point lies close to the neutron star. 

We note that the discussion in this section is valid only if the magnetic fields are closed 
throughout the disk. Only in these cases does Iq measure the net torque on the neutron star. Thus 
our results in this section are much more restrictive compared to the location of the sonic point 
(which depends only on the local field structure) discussed in §§3-4. 



6. Discussion: Applications to QPOs in LMXBs 



We have presented in §§2-4 an unified treatment of the inner region of accretion flow under 
the combined influences of general relativistic gravity and stellar magnetic fields in LMXBs. We 
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have shown that even relatively weak magnetic fields (10^ — 10^ G) can slow down the orbital 
motion in the inner disk by taking away angular momentum from the disk, thereby changing the 
position of the sonic point significantly (cf. Figs. 2-3). 

While the mechanisms responsible for the kHz QPOs in LMXBs arc still uncertain, it is 
tempting to associate them with orbital motions at a certain preferred radius. If this is the case, 
then the Keplerian frequency at inner radius of the disk, or more precisely the sonic radius is 
certainly a natural choice (see Paczynski 1987 for an earlier suggestion on the importance of the 
disk sonic point). One can envisage a number of different mechanisms that will lead to QPOs at 
uxirs) — i^so and the beat frequency with the stellar spin (Miller et al. 1996; Van der Klis 1997). 
The following discussion is based on the hypothesis that the kHz QPO frequency corresponds to 
i^MSO or the sonic-point Kepler frequency. We note that although our treatment of the magnetic 
field effects in §2-4 is rather general (albeit based on a phenomenological prescription), other 
physical effects might be at work (such as radiation forces for high-luminosity systems; see §1). As 
a result, some of our conclusions below should be considered tentative. 

(a) Range of QPO frequencies and constraint on the magnetic field strength: As emphasized by 
Van der Klis (1997), similar QPO frequencies (500 — 1200 Hz) have been observed in sources with 
widely different average luminosities {L) (from a few times 10""^LEdd to near L^dd, corresponding 
to (M) between a few times 10^^ g s~^ to 10^^ g s~^), while for an individual source i^qpo 
often correlates strongly with the X-ray flux. This peculiar lack of correlation between the QPO 
frequency and (L) can be explained in our model, as long as the star's magnetic field strength 
correlates with its (M) in such a way as to leave {M) /{(jB^) in a certain range (see Fig. 3). 
For example, if we use the model with n = 3 (dipole field) and Vg = 300 Hz, then to produce 
UMSO in the range of 500 - 1200 Hz would require 0.003 <, Mi7/(/3S|) <, 0.1; on the other 
hand, if we use the model with n = 2 and 1/5 = (open field configuration), we would require 
0.05 ^ Mi7/{(3B^) ^ 0.3 (we have adopted Mi. 4 = Riq = 1 in these examples). Indeed, such 
correlation between the magnetic field strength and the mean mass accretion rate among LMXBs 
has been suggested independently on the basis of Z and atoll source phenomenology (Hasinger &; 
Van der Klis 1989), although the origin for this correlation is still unclear. We note, however, 
that the correlation needs not be very strong, considering the wide range of other controlling 
parameters such as n (the shape of the magnetic field) and Ug (stellar rotation) (see Fig. 3). 

While there is a natural upper limit to z^so (corresponding to tmso 6M; but see (c) 
below), the existence of a (source-dependent) lower limit to the observed QPO frequency needs an 
explanation. If we rely on inner disk accretion to explain these QPOs (e.g.. Van der Klis 1997), 
then one possibility is that for large tmso (small z^mso); the accreting gas can be channeled out of 
the disk plane by the magnetic field — This must happen for sufficiently small M (or sufficiently 
large 5, as in the case of accreting X-ray pulsars). The precise location where the plasma leaves 
the disk depends on the near-zone field structure, and is clearly source-dependent. 

The observed QPO frequencies may already be used to probe the magnetic fields in LMXBs 
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and the nature of the magnetic field - disk interactions. As an example, consider the atoll source 
4U 0614+091 which has extremely small luminosity (see Ford et al. 1997 and references therein): 
To obtain i'mso ~ 1150 Hz requires Mn/{(3Bj) ^ 0.1 (This constraint depends somewhat on the 
field structure and Vg', see Fig. 3). At M17 ~ 0.05, this translates to P^^'^Bj 0.7 — a stronger 
magnetic field would push the sonic point (or the generalized MSO) to a larger radius. We are left 
with two possibilities: (i) If Bj ^ 10, we would require 13 <^ 0.01, i.e., the dissipation of toroidal 
fields near the disk must be very efficient; (ii) If /3 ~ 1 (as expected for open field configurations), 
then we would require Bj ^ 1. Indeed, if the kHz QPO sources represent a fair sample of LMXBs, 
then we might conclude that the magnetic fields in LMXBs are systematically weaker than those 
in millisecond pulsars. This may indicate that the magnetic field of a neutron star is "buried" 
during the LMXB phase (e.g., Romani 1990; Urpin & Geppert 1995; Konar &: Bhattacharya 1997; 
Brown & Bildsten 1998), and later regenerates or re-emerges as accretion stops. 

(b) Scaling of v'qpo with M: For most sources^, it was found that the kHz QPO frequency 
strongly correlates with the XTE count rate (2 — 50 keV), with power-law index greater than 
unity. However, the scaling relation between the count rate and M is not well established, and it 
has been suggested the flux of the black-body component is a better indicator of QPO frequency 
(Ford et al. 1997). As discussed in §3, i^mso depends primarily on the magnetic field structure 
near the sonic point, particularly on the "field shape" index n. If the scaling of z^qpo with M can 
be established observationally, it may be possible to distinguish a closed field configuration from 
an open one. For example, if we believe the scaling z^qpo oc M'^ with a > 1, then we may conclude 
that n < 2 (see Eq. and the discussion following Eq. [^]), which indicates that magnetic 
fields in LMXBs do not have dipolar shape, but rather have complex topology (see Arons 1993 for 
discussion on related issues). 

(c) The maximum value of vqpq: It has been suggested (Zhang et al. 1997) based on the 
narrow range of the maximal QPO frequencies (1100 — 1200 Hz) in at least six sources that these 
maximum frequencies correspond to the Kepler frequency at tgr = 6M, which then implies 
that the neutron star masses are near 2Mq (see also Kaaret et al. 1997). While we agree that 
this conclusion seems most natural, we nevertheless add the following cautionary notes: (i) The 
inferred large stellar masses may be problematic: All neutron stars with well-determined masses 
(including a few that certainly had accreted mass, although not necessarily in the same accretion 
mode as in LMXBs) have masses consistent with being M ~ 1.4 (e.g., van Kerkwijk et 

al. 1995). In particular, the 5.4 ms recycled pulsar B1855+09, which is thought to have gone 
through a LMXB phase (Phinney & Kulkarni 1994), has a mass 1.50 ±[];f^ M© (Kaspi et al. 1994). 
Moreover, accretion of 0.6Mq might have spun up the neutron stars to near break-up (see Cook 
et al. 1994 for calculations of spin- up tracks in the nonmagnetic case), in contrary to the observed 



*The possible exceptions are 4U 1608-52 (Berger et al. 1997) and the high-intensity ("banana") state of 4U 1636-53 
(Wijnands et al. 1997), for which the power spectrum shows a single peak with frequency independent of the count 
rate. However, this may be an artifact of insensitive search of weak QPOs; see Van der Klis (1997). 
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spin rates (300 — 500 Hz), (ii) If the maximum vqpq is indeed i/i<:(6M), then the correlation 

between fqpo and M should weaken as M increases, and eventually i^qpo should approaches a 
constant independent of M (see Fig. 3). This has not been observed. Therefore in our opinion it is 
premature to identify the maximal fQpo with the Kepler frequency at 6M. An alternative is that 
as ^Qpo approaches 1100 — 1200 Hz, the rms QPO amplitude decreases — as the observations 
have indicated, making it difficult to detect higher QPO frequencies. 

(d) Horizontal- Branch Oscillations (HBOs) in Z-sources: In several Z-sources (e.g., Sco X-1, 
GX 5-1 and GX 17+2), HBOs with frequencies 20 — 50 Hz have been detected simultaneously 
with the kHz QPOs (see Van der Klis 1997). One standard interpretation of HBOs is that they 
are associated with the beat between the Kepler frequency at the magnetosphere boundary and 
the neutron star spin (Alpar & Shaham 1985). Since the spin frequencies i^g of these sources (as 
determined from the difference in the twin kHz QPO frequencies) lie around 300 Hz (see White 
Sz Zhang 1997), the putative magnetosphere boundary must be located at a large radius where 
i^K ^ i^s ^ 300 Hz. As we have shown in this paper, such a strong magnetic field must necessarily 
push the (generalized) disk sonic point to a large radius where the Kepler frequency drops below 
the kilo-Hertz range. (Recall that for low field systems such as LMXBs, the distinction between 
the sonic point and the magnetosphere boundary probably does not exist, and the two separate 
radii are replaced by a single generalized sonic point [see §§3-4].) Therefore if the kHz QPOs are 
associated with the sonic-point Kepler frequency, then the magnetospheric beat frequency model 
for HBOs cannot work, and the origin of HBOs must lie elsewhere. Alternative models for HBOs 
have been discussed by Biehle & Blandford (1993) and Stella & Vietri (1997). 

The author thanks Zhiyun Li, Rob Nelson and Brian Vaughan for valuable discussions and 
Lars Bildsten for comment. He also thanks the referee for constructive comments which improved 
the presentation of the paper. This research is supported by a Richard Chace Tolman Fellowship 
at Caltech, NASA Grant NAG 5-2756, and NSF Grant AST-9417371. 
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Fig. 1. — (a) Radial velocity \u\ (in units of c) and (b) angular velocity Q, (in units of c^/GM) 
of transonic accretion flows as a function of disk radii. The model parameters are n = 3 (dipolar 
field), R/M = 5, a = 0.1, Cg = 0.01, = 0.013/M, and 6^ = i (solid curves), 10 (dashed curves). 
The sonic points are marked by filled circles. The dotted lines depict the sound speed (a) and 
Keplerian angular velocity (b). 
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Fig. 2. — (a) The sonic radius and (b) the constant angular momentum as a function of the 
dimensionless ratio 6^ (defined in Eq. [18]) for three different models: The dotted lines are for 
a = 0.1, Cs = 0.01, the short-dashed lines for a = 0.02, Cg = 0.01, and the long-dashed lines for 
a = 0.02, Cs = 0.005. All models have n = 3, R/M = 5 and Vis = 0.013/M. The sohd lines are the 
asymptotic analytic solutions as discussed in §4. 
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Fig. 3. — Orbital frequency z^so = I'Ki'r'MSo) at the generalized MSO (which approximates the 
sonic point) as a function of mass accretion rate. The results are obtained using the analytical 
expressions given in §4. The dotted line corresponds to an open field configuration with fig = 
and n = 2, while the other three lines correspond to closed field configurations with n = 3, and 
fls = (solid line), 0.013/M (dashed line), and 0.026/M (long-dashed line). Note that for large 
M/{PB^) (or small b^), ruso approaches tgr = 6M and z^mso approaches 1.57/Mi.4 kHz. 
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Fig. 4. — The equilibrium spin frequency i^g of the neutron star as a function of the dimensionless 
ratio 6^ (defined in Eq. [18]). The sohd Une is our result including both the GR and magnetic field 
effects, and the dashed line is usual result as given by Eq. (36). 



